A new measurement of the altitude dependence of the atmospheric muon intensity 
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We present a new measurement of atmospheric muons made during an ascent of the High Energy 
Antimatter Telescope balloon experiment. The muon charge ratio /x + j pT as a function of atmo- 
spheric depth in the momentum interval 0.3-0.9 GeV/c is presented. The differential pT intensities 
in the 0.3-50 GeV/c range and for atmospheric depths between 4-960 g/cm 2 are also presented. 
We compare these results with other measurements and model predictions. We find that our charge 
ratio is ~1.1 for all atmospheric depths and is consistent, within errors, with other measurements 
and the model predictions. We find that our measured \i~ intensities are also consistent with other 
measurements, and with the model predictions, except at shallow atmospheric depths. 
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I. INTRODUCTION 



Measurements and theoretical calculations of atmo- 
spheric neutrino fluxes have consistently disagreed, and 
this disagreement was interpreted in terms of neutrino os- 
cillations. However, this interpretation requires an accu- 
rate understanding of neutrino fluxes in the atmosphere. 
Simulations of the absolute fluxes have been conducted 
(JUS El and references therein) but suffer from system- 
atic uncertainty in the normalization of the atmospheric 
neutrino spectrum. These simulations also predict the 
spectrum of other atmospheric secondaries, specifically 
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muons. Indeed, the production of neutrinos in the atmo- 
sphere is closely coupled with that of muons, as they are 
produced together in pion and kaon decays, and as some 
muons themselves decay and contribute to the neutrino 
flux. Therefore, measurements of atmospheric muons by 
ascending high altitude balloon borne instruments can 
be used to reduce the neutrino model uncertainties. It 
should be pointed out that other experimentally accessi- 
ble quantities can also be used to that end, such as the 
absolute primary cosmic ray flux impinging on the atmo- 
sphere, or other atmospheric secondaries such as gamma 
rays or secondary nucleons. The atmospheric muon flux 
is therefore an important element in any comprehensive 
strategy to improve our understanding of neutrino pro- 
duction in the atmosphere. 

An earlier version of the High Energy Antimatter Tele- 
scope (HEAT-e^ instrument had been used to measure 
air-shower muons during atmospheric ascent \4. The 
HEAT instrument is described elsewhere 0, 07 In its 
present configuration, HEAT-pbar is optimized to study 
antiprotons. It combines a superconducting magnet spec- 
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trometer using a drift-tube hodoscope (DTH), a time-of- 
fiight system (TOF), and two stacks of multiwire propor- 
tional chambers (dE/dx). We report here a new mea- 
surement of the muon charge ratio as a func- 
tion of atmospheric depth in the momentum interval 0.3- 
0.9 GeV/c, and differential fi~ intensities in the 0.3-50 
GeV/c range and for atmospheric depths between 4-960 
g/cm 2 for a balloon flight from Fort Sumner, NM, USA 
on June 3, 2000. 



II. MUON IDENTIFICATION AND 
BACKGROUNDS 

The ascent phase of the flight lasted ~3 hours during 
which more than 1.5 million events were recorded at an 
average vertical geomagnetic rigidity cutoff of 4.5 GV. 
The flight occurred under solar maximum conditions (<!> 
= 1330 MV). The atmospheric overburden ranged from 
a high of 960 g/cm 2 with the instrument on the ground 
to a low of 4 g/cm 2 at float altitude. 

Muon events were identified using the three main de- 
tectors onboard the HEAT-pbar instrument. The meth- 
ods of muon identification and selection criteria involving 
the TOF and DTH detectors are described in detail else- 
where |4|. To select muon events using the TOF detector, 
the incident particle must have 0.85 < (3 < 2.00, where 
(3 — v/c is the particle's velocity in units of the speed of 
light. Also, for the particle must have 0.3 GV < R 
< 0.9 GV and for fi~, R < -0.3 GV, where R = pc/Ze 
is the particle's magnetic rigidity, or momentum per unit 
electric charge. Muon events with R > 0.9 GV cannot be 
selected based on rigidity as a function of /3, due to pro- 
ton contamination. The range of selected also ensures 
that the particle was traveling downward through the in- 
strument. To select muon events with the DTH detector, 
the particle track must have passed through at least 4 of 
the 8 drift tube planes in the non-bending view, and at 
least 8 of the 18 planes in the bending view. This is to 
ensure that there are enough samples along the particle's 
track to produce a valid track and rigidity measurement 
through the DTH. We retain only events with average 
fit residuals along the track that are better than 6 mm 
in the x direction, 1.5 mm in the y direction, and 1 mm 
in the z direction. This rejects only ^0.1% of events, 
and provides basic track cleanliness criteria. Lastly, the 
event is required to have a tracking fit goodness parame- 
ter Xdt — 1^.0 and | (Maximum Detectable Rigidity)/R| 
> 4, which are used to ensure good track quality. 



A. The dE/dx Detector 

1. Physical Construction 

The dE/dx detector consists of two stacks of 70 seg- 
mented multi-wire proportional chambers, one stack 
above and one below the magnetic hodoscope as shown 



in Figure^ The 140 chambers are each ~1 cm thick and 
filled with a gas mixture of 95% Xenon and 5% Methane 
(used as a quench gas). Xenon was used since it yields 
the strongest possible logarithmic rise in the energy loss 
curve from minimum ionization to relativistic saturation 
(70%), and hence gives the highest possible velocity res- 
olution. 




FIG. 1: Schematic view of the HEAT-pbar Instrument. The 
solid black line labeled "S" is a piece of scintillator and is part 
of the instrument trigger. 



2. Energy Loss Determination 

Each chamber's anode wires are connected in groups of 
4-6 wires. In addition to the ionization energy deposits, 
there are occasional spurious signals in wire groups away 
from the particle track. For analysis, only wire groups 
less than 5 cm (the average width of two adjacent wire 
groups) from the particle's track are used in the energy 
loss determination for each chamber. The signals are 
corrected for detector response changes during the flight, 
electronic gain variations, and particle incident angle. 
Some particle trajectories are within the instrumental 
acceptance but miss some of the chambers near their 
edges, so that the mean number of chambers hit in an 
event is 128 out of 140. We require the particle to have 
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left a signal in at least 103 of 140 chambers to ensure 
sufficient sampling points for a determination of the av- 
erage energy loss. The restricted mean average energy 
loss is determined by eliminating chambers whose signal 
is not in the lowest 58% of the signals from all chambers 
hit. The restricted mean average is used to minimize the 
effects of Landau fluctuations on the final energy loss de- 
termination. A typical plot of restricted mean energy loss 
(dE/dx) versus positive rigidity for ascent data is shown 
in Figure [21 Included on the plot are the Bethe-Bloch 
energy loss curves for deuterons, protons, muon/pions, 
and positrons. Various particle species clearly populate 
the plot. Note the geomagnetic cutoff effect on the pro- 
ton data at -4.5 GV (Fort Sumner, NM, USA). FigureEl 




Log(Rigidity) in GV 



FIG. 2: Distribution of dE/dx as a function of rigidity for as- 
cent data for positive rigidities. Superimposed are the Bethe- 
Bloch energy loss curves for various particle species. 

shows the restricted mean energy loss for events in the 
rigidity range -0.9 GV < R < -0.3 GV, as an example, 
which highlights the highly gaussian nature of the re- 
stricted mean energy loss. Distributions in other rigidity 
ranges look similar. 



3. dE/dx Event Selection 
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Restricted Mean Energy Loss (dE/dx) in arbitrary units 

FIG. 3: Restricted mean energy loss distributions for ascent 
data with gaussian fits for -0.9 GV < R < -0.3 GV. 



restricted mean energy loss versus rigidity. 
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For the dE/dx detector, the difference between the re- 
stricted mean energy loss for the upper and lower dE/dx 
stacks must be within approximately two standard devi- 
ations from the mean of the difference distribution. This 
ensures that an event that showers below the top dE/dx 
stack is rejected. The restricted mean dE/dx signal of 
the particle must fall between and 65 (inclusive, in ar- 
bitrary units) for 0.3 GV < |R| < 0.8 GV, between and 
70 for 0.8 GV < |R| < 4.0 GV, and between and 78 
for 4.0 GV < |R| < 50.0 GV. See Figure H for a plot of 



FIG. 4: Restricted mean energy loss as a function of rigidity 
for ascent data for (top) positive rigidities and (bottom) neg- 
ative rigidities. The event selections (solid lines) are shown 
on each plot as well as the Bethe-Bloch energy loss predic- 
tion curves (dashed). Note that the apparent higher density 
of negative muons compared to positive muons is an artifact 
of the plotting routine used for the figure; there are approxi- 
mately equal numbers of points shown on each panel. 
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B. Backgrounds to the Differential Muon Intensity 

Protons, electrons, and pions are the major back- 
ground contaminations of the differential muon intensity. 
The pion intensity, though relatively small with respect 
to the muon intensity, increases with altitude. From 0, 
the pion/muon ratio is —0.035 at a depth of 3 g/cm 2 . 
Fortunately, the intensities of other products of atmo- 
spheric showers, such as kaons, are small in relation to 
the muon intensity. 

Particle interactions within the HEAT-pbar instru- 
ment also produce pions that can contaminate the muon 
intensity. A GEANT 3.21 Monte Carlo simulation of the 
instrument with an input proton spectrum falling with an 
index of —2.7 and an input momentum range of 4.5 GV < 
R < 50 GV yielded a secondary interaction pion to pro- 
ton fraction of 1.1 x 1CT 3 , which decreased to 9.7 x 1CT 5 
after application of the muon selection criteria. In or- 
der to determine what percentage of the muon intensity 
this represents, and hence the contamination percentage, 
the proton to muon ratio must be known. As a first or- 
der estimate, we used proton and muon information from 
Q , and determined that the proton-to-muon ratio ranges 
from 1.6 x 10~ 4 at the ground to a maximum of 3.0 at 
13 g/cm 2 . Choosing the highest value for the ratio of- 
fers a worst case contamination fraction, since the ratio 
decreases with depth. Using a ratio of 3.0 yielded an 
interaction pion intensity that was 0.01% of the proton 
intensity and 0.03% of the muon intensity at atmospheric 
depths greater than 13 g/cm 2 . These percentages only 
decrease further with atmospheric depth, so the interac- 
tion pion contamination of the muon intensity is ignored 
since the error introduced is much less than the statis- 
tical error of the measurement. The number of kaons 
produced in the instrument is even smaller than the pion 
number, so they are ignored as well. 

III. RESULTS AND DISCUSSION 

A. Muon Charge Ratio 

The number of muons for various atmospheric depths 
for the rigidity range 0.3 < |R| < 0.9 GV is shown in 
Table 0] The /i + / fi~ ratio is shown as a function of 
atmospheric depth in Figure [S] for a geomagnetic cut- 
off of -4.5 GV (Fort Sumner, NM, USA), and in Fig- 
ure [fjj for a very small cutoff at northern latitudes (Lynn 
Lake, Manitoba, Canada). Shown are the new HEAT- 
pbar measurements (labeled HEAT 00), as well as two 
previous HEAT-e measurements and other recent mea- 
surements The various measurements were 
made at different solar epochs, so that direct compari- 
son is problematic. However, there is good agreement 
between the various measurements at Fort Sumner, irre- 
spective of the solar modulation conditions, which is to 
be expected since the primary flux modulation is small 
above 4.5 GV. There are only two measurements at Lynn 



TABLE I: Muon charge ratio as a function of atmospheric 
depth (d) for 0.3 < |R| < 0.9 GV. Errors on the ratios are 
statistical only. 



d 


d 


N + 


N - 

M 


N +/N - 


(g/cm 2 ) 


(g/cm 2 ) 








4-7 


5.572 


228 


200 


1.14 ± 0.11 


7-13 


9.583 


428 


346 


1.237 ± 0.089 


13-32 


22.49 


910 


837 


1.087 ± 0.052 


32-67 


49.68 


1539 


1251 


1.230 ± 0.047 


67-140 


105.2 


2460 


2099 


1.172 ± 0.035 


140-250 


190.8 


1552 


1447 


1.073 ± 0.039 


250-350 


296.7 


658 


642 


1.025 ± 0.057 


350-840 


497.5 


949 


852 


1.114 ± 0.053 


840-960 


887.3 


688 


562 


1.224 ± 0.070 



Lake, with relatively large uncertainties, and where solar 
modulation effects are expected to be more important. 
Both measurements were made under near solar max- 
imum conditions, and are compatible with each other 
given the large uncertainties. Also shown on Figures [5] 
and are 3-dimensional calculations of the air shower 
development with the TARGET algorithm [n), widely 
used for neutrino flux calculations. Curves are shown for 
both solar maximum and minimum conditions, for aver- 
age primary fluxes (which may not represent the actual 
spectrum at the time of any individual measurement). 
The agreement between the calculations and the mea- 
surements (especially the HEAT measurements) appears 
reasonably good, with the caveat that statistical uncer- 
tainties at high altitudes for the HEAT 95 measurements 
are large. 

B. /i~ Differential Intensity 

Our measurement of the fj,~ differential intensities was 
made using similar corrections as for previous BEAT- 
ON flights with only minor modifications for differing 
instrument configuration. The results are shown in Fig- 
ures □ HI and El and listed in Table HU The differential 
intensities as a function of rigidity and atmospheric depth 
are shown multiplied by R 2 for better comparison. Also 
shown are the results of other measurements f|9l ll0lll2| ). 
The level of agreement is quite good despite the differing 
solar epochs and rigidity cutoffs. The curves shown are 
for calculations with the TARGET algorithm with the 
same parameters as for Figure |3J 

The agreement between the HEAT measurements and 
model predictions is for the most part excellent, except 
for the very highest altitudes where the predictions are 
below the measured intensities. We note that in Fig. [7| 
there appears to be a tendency for the measurements of 
the ground /i~ intensities made at a 4.5 GV geomagnetic 
cutoff to be slightly higher than those at sub-GV cutoff. 
We do not expect this effect to be real, as energetic muons 
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FIG. 5: Atmospheric muon charge ratio as a function of at- 
mospheric depth for a geomagnetic cutoff of 4.5 GV (Fort 
Sumner, NM, USA). Also shown are the ratios for the HEAT 
1994 flight as well as other measurements. The curves are 
calculations with the 3-dimensional TARGET algorithm. In- 
dicated are the solar modulation parameters for each data 
set. 



(R > 1 GV) at the ground are the end products of parent 
cosmic rays well above the geomagnetic cutoff at either 
location. We note that differential intensities are rescaled 
by R 2 , so that the effect of any small systematic error in 
the rigidity determination is magnified. 



IV. CONCLUSIONS 

We have measured with good statistics the muon 
charge ratio /i + / [i~ as a function of atmospheric depth 
in the momentum interval 0.3-0.9 GeV/c and the dif- 
ferential /i~ intensities in the 0.3-50 GeV/c range and 
for atmospheric depths between 4-960 g/cm 2 . We have 
found that our charge ratio is ~1.1 for all atmospheric 
depths and is consistent, within errors, with other mea- 
surements and the model predictions. We have found 
that our measured jjT intensities are also largely consis- 
tent with other measurements, and with the model pre- 
dictions, despite varying solar epochs and geomagnetic 
rigidity cutoffs. This reinforces the conclusion of 4] that 
model calculations are adequate to the task of predicting 
neutrino fluxes. The discrepancy between our measure- 
ments and the model predictions at shallow atmospheric 
depths is of little import to our understanding of neutrino 
fluxes, as the bulk of the neutrinos come from decay pro- 
cesses much deeper in the atmosphere. 



FIG. 6: Atmospheric muon charge ratio as a function of at- 
mospheric depth for a geomagnetic cutoff of <C 1 GV (Lynn 
Lake, Manitoba, Canada) for two data sets. The curves are 
calculations with the 3-dimensional TARGET algorithm. In- 
dicated are the solar modulation parameters for each data 
set. 




1 10 

Rigidity (GV) 

FIG. 7: Ground level differential fi~ intensities as a function 
of rigidity. Also shown are results from various other instru- 
ments. The plots have been scaled by a factor of R 2 . The 
curves are calculations with the 3-dimensional TARGET algo- 
rithm for solar minimum (solid) and solar maximum (dashed). 
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FIG. 8: Differential /i~ intensities as a function of rigidity for various atmospheric depth ranges. Also shown are results 
from various other instruments. (Note: only subsets of measurements by other groups are shown here, where the choice of 
atmospheric depth bins matches with ours.) The plots have been scaled by a factor of R 2 . The curves are calculations with 
the 3-dimensional TARGET algorithm for solar minimum (solid) and solar maximum (dashed). See Figures |S] and |5] for each 
measurement's solar and geomagnetic conditions. 
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FIG. 9: fi~ growth curves for various rigidity intervals. (Note: other groups have used different rigidity bins in their analyses, 
so that we only report HEAT measurements here.) The plots have been scaled by a factor of R 2 . The curves are calculations 
with the 3-dimensional TARGET algorithm for solar minimum (solid) and solar maximum (dashed). See Figure Q for each 
measurement's solar and geomagnetic conditions. 
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TABLE II: Atmospheric fj, intensities [in (cm 2 s sr GeV/c) x ] and number of recorded events as a function of rigidity and 
atmospheric depth. 





4-7 g/cm 2 


7-13 g/cm 2 




13-32 g/cm 2 


32-67 g/cm 2 


67-140 g/cm 2 




0.3-0.5 


(1.51 ± 0.26) x 10~ 3 


(2.06 ± 0.31) x 10 


-3 


(4.62 ± 0.61) x 10~ 3 


(8.5 ±1.1) x 10 -3 


(1.18 ± 0.15) x 10" 


-2 


GV 


61 


108 




279 


415 


672 




0.5-0.8 


(1.23 ± 0.18) x 10~ 3 


(1.58 ± 0.21) x 10 


-3 


(3.27 ± 0.41) x 10~ 3 


(6.14 ± 0.74) x 10 -3 


(8.82 ± 1.09) x 10 


-3 


GV 


110 


186 




440 


670 


1116 




0.8-1.2 


(5.21 ± 0.85) x 10~ 4 


(1.02 ± 0.14) x 10 


-3 


(1.84 ± 0.24) x 10~ 3 


(3.61 ± 0.44) x 10~ 3 


(5.58 ± 0.69) x 10 


-3 


GV 


71 


185 




380 


607 


1086 




1.2-1.7 


(4.08 ± 0.67) x 10~ 4 


(5.26 ± 0.76) x 10 


-4 


(1.12 ± 0.15) x 10~ 3 


(2.32 ± 0.29) x 10~ 3 


(3.10 ± 0.39) x 10 


-3 


GV 


71 


120 




294 


498 


768 




1.7-2.5 


(-i 40 + n 27) x 10 -4 


(o 7K -1- n 41 \ x 1(1 


-4 


(c. co 1 n x in-4 


( 1 ia 4- n 1 K\ x 10 — 3 




-3 


GV 


41 


104 




246 


415 


629 




2.5-4.0 


(6.7 ± 1.3) x 10" 5 


(1.01 ±0.16) x 10 


-4 


(2.10 ±0.29) x 10~ 4 


(4.63 ±0.59) x 10~ 4 


(6.34 ±0.81) x 10 


-4 


GV 


37 


73 




174 


311 


494 




4.0-8.0 


(3.11 ±0.59) x 10~ 5 


(4.13 ±0.67) x 10 


—5 


(8.0 ± 1.1) x 10" s 


(1.43 ±0.19) x 10~ 4 


(1.88 ±0.24) x 10 


-4 


GV 


43 


75 




165 


239 


363 




8.0-16.0 


(6.6 ± 1.8) x 10~ 6 


(9.9 ±2.0) x 10" 


6 


(1.62 ±0.28) x 10" 6 


(2.87 ± 0.45) x 10" 6 


(4.01 ±0.58) x 10 


-5 


GV 


17 


34 




63 


90 


145 




16.0-32.0 


(1.52 ±0.64) x 10~ 6 


(1.68 ±0.62) x 10 


-6 


(4.3 ± 1.0) x lO" 6 


(6.1 ± 1.4) x lO -6 


(7.2 ± 1.4) x 10~ 


6 


GV 


6 


8 




23 


25 


39 




32.0-50.0 




(4.1 ±4.1) x 10" 


7 


(3.7 ±3.7) x 10~ 7 


(5.4 ±5.4) x 10~ 7 


(3.0 ± 1.1) x 10~ 


6 


GV 





1 




1 


1 
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140-250 g/cm 2 




250-350 g/cm 2 


350-840 g/cm 2 




840-960 g/cm 2 




0.3-0.5 


(1.88 ±0.24) x 10" 


-2 


(1.40 ±0.19) x 10- 2 


(6.69 ± 0.88) x 10" 


-3 


(2.28 ±0.31) x 10- 


-3 


GV 


477 




213 


247 




172 




0.5-0.8 


(1.29 ±0.16) x 10 


-2 


(9.7 ± 1.2) x 10- 3 


(5.51 ± 0.68) x 10" 


-3 


(1.64 ±0.21) x 10" 


-3 


GV 


730 




329 


457 




283 




0.8-1.2 


(8.4 ± 1.1) x 10- 


3 


(6.83 ± 0.86) x 10" 3 


(3.97 ±0.49) x 10" 


-3 


(1.39 ±0.18) x 10" 


-3 


GV 


737 




359 


509 




373 




1.2-1.7 


(4.84 ±0.62) x 10" 


-3 


(4.37 ±0.56) x 10- 3 


(2.74 ±0.34) x 10" 


-3 


(1.22 ±0.15) x 10" 


-3 


GV 


538 




291 


448 




421 




1.7-2.5 


(2.72 ±0.35) x 10 


-3 


(2.23 ±0.29) x 10~ 3 


(1.47 ±0.18) x 10" 


-3 


(8.01 ±0.99) x 10" 


-4 


GV 


501 




247 


396 




457 




2.5-4.0 


(1.00 ±0.13) x 10 


-3 


(8.2 ± 1.1) x 10- 4 


(7.57 ±0.95) x 10- 


-4 


(4.37 ± 0.54) x 10" 


-4 


GV 


348 




171 


392 




482 




4.0-8.0 


(3.34 ±0.45) x 10" 


-4 


(2.67 ±0.38) x 10- 4 


(2.69 ± 0.34) x 10~ 


-4 


(1.77 ±0.22) x 10- 


-4 


GV 


288 




138 


349 




491 




8.0-16.0 


(5.89 ±0.93) x 10 


-5 


(5.27 ± 0.95) x 10" 5 


(6.20 ±0.87) x 10- 


-5 


(4.64 ± 0.62) x 10- 


-5 


GV 


95 




51 


152 




230 




16.0-32.0 


(9.5 ±2.4) x 10" 


6 


(8.4 ±2.6) x 10" 6 


(9.0 ± 1.9) x 10-' 




(1.15 ±0.18) x 10- 


-5 


GV 


21 




12 


30 




93 




32.0-50.0 


(1.9 ± 1.4) x 10- 


6 


(3.9 ± 2.3) x 10- 6 


(3.0 ± 1.4) x 10-' 




(2.66 ±0.83) x 10- 


-6 


GV 


2 




3 


5 




12 





